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We use the non-equilibrium Green's function formalism and a self-consistent Hartree-Fock approx- 
imation to numerically study the effects of a single impurity and interactions between the electrons 
(with and without spin) on the conductance of a quantum wire. We study how the conductance 
varies with the wire length, the temperature, and the strengths of the impurity and interactions. 
The numerical results for the dependence of the conductance on the wire length and temperature 
are compared with the results obtained from a renormalization group analysis based on the Hartree- 
Fock approximation. For the spin-1/2 model with a repulsive on-site interaction or the spinless 
model with an attractive nearest neighbor interaction, we find that the conductance increases with 
increasing wire length or decreasing temperature. This can be explained using the Born approx- 
imation in scattering theory. For a strong impurity, the conductance is significantly different for 
a repulsive and an attractive impurity; this is due to the existence of a bound state in the latter 
case. In general, the large density deviations close to the impurity have an appreciable effect on the 
conductance at short distances which is not captured by the renormalization group equations. 



PACS number: 73.23.-b, 73.63. Nm, 71.10.Pm 



I. INTRODUCTION 

The conductance of electrons in a quantum wire has 
been the subject of intensive study in recent years, both 
experimentally [1-6] and theoretically [7-9]. For a wire in 
which only one channel is available to the electrons and 
the transport is ballistic (i.e., there are no impurities in- 
side the wire, and there is no scattering from phonons or 
from the contacts between the wire and its two leads), 
the conductance is given by G — 2e 2 /h for infinitesi- 
mal bias [10,11]. This result is expected to hold even 
if one takes into account the interactions between the 
electrons since such two-body scatterings conserve the 
momentum. However, if there is an impurity inside the 
wire which scatters the electrons, then the conductance 
is reduced because such a scattering does not conserve 
the momentum of the electron. For a one-dimensional 
system containing a ^-function impurity with strength 
V, we obtain 

G = ^(1 - cV 2 ), (1) 

to lowest order in V, where c is a constant related to the 
Fermi velocity of the electrons (see the discussion below 
Eq. (3)). (The situation is different if the wire is only 
quasi-one-dimensional, and the impurity potential has a 
finite range [12]). In the absence of interactions, G does 
not depend on the wire length L or the temperature T 
(as long as ksT is much less than the Fermi energy). 
But in the presence of interactions, it turns out that V 
effectively becomes a function of the length scale (which 
is related to either L or T as will be explained below), 
and G therefore varies with L and T. The variation of V 



with length scale is governed by a renormalization group 
(RG) equation. (In this paper, we will only consider a 
non-magnetic impurity which scatters electrons in a spin 
independent way). 

There are three length scales which are of interest in 
the problem. The smallest of them is A = ir/kp (where 
kp is the Fermi wavenumber); this is the wavelength of 
the oscillations in the electronic density near impurities 
as we will see. The other two length scales are the length 
of the wire L, and the thermal coherence length Lp which 
is equal to Tivf / '(ksT) (where vf is the Fermi velocity). 
Lt gives an idea of the distance beyond which an electron 
wave function loses its phase coherence. At a tempera- 
ture T, the conductance typically receives contributions 
from a number of states near the Fermi energy whose en- 
ergies have a spread of the order of AE = ksT. Hence 
the spread in momentum Ap = AE/vf is of the order 
of ksT/vp- A superposition of waves with such a spread 
of momenta loses phase coherence in a distance of the 
order of Lp = %vf / (kpjT). Electronic transport is there- 
fore thermally incoherent if the wire length L >> hp, 
coherent if L << Lp, and partially coherent in the inter- 
mediate range. This will become clearer when we discuss 
our numerical results. 

The RG equations for the transmission coefficient \t\ 2 
has been derived from continuum theories in several ways 
[7,8,13-15]. We will discuss these equations for two mod- 
els in Sec. II. (The conductance is related to the trans- 
mission as G = (2e 2 / h)\t\ 2 , where the factor of 2 is due 
to the electron spin). The RG equation is used to ana- 
lytically follow the evolution of |<| 2 starting from a short 
distance scale and going up to a length scale which is the 
smaller of the two quantities L and Lp. If L >> Lp, 



the conductance is governed by Lt and not L, and vice 
versa if Lt >> L. The RG equation has two fixed points 
which lie at \t\ 2 = and 1; the system approaches one of 
these fixed points if L and Lp are both very large. 

The above statements are only valid for length scales 
much longer than A because only then can one use a 
continuum description from which the RG equation is de- 
rived. It would therefore be useful to consider an alterna- 
tive method for computing the conductance which works 
even at short length scales where the continuum descrip- 
tion is not valid. An interesting thing which occurs at 
short distances is that if the impurity provides an at- 
tractive potential to the electrons, there is a bound state 
whose wave function decays exponentially away from the 
impurity. We may wonder what effect a bound state has 
on the conductance; the RG equations mentioned above 
do not take this state into account. One may think that 
a bound state (whose energy lies outside the bandwidth 
of the leads) cannot directly affect the conductance, be- 
cause electrons coming in from or going out to the leads 
cannot enter or leave such states in the absence of any 
inelastic scattering. However, such a state contributes 
to the electronic density near the impurity, and that can 
affect the transport in the presence of interactions. 

In a series of papers, the technique of functional RG 
has been used to numerically study the conductance and 
other properties of interacting electron systems in one di- 
mension [16-19]. The authors of those papers have gone 
up to very large system sizes and have found excellent 
agreement at those length scales between their numeri- 
cal results and the asymptotic scaling forms given by the 
RG equations. 

In this paper, we will use the non-equilibrium Green's 
function (NEGF) formalism to numerically study the 
conductance of a quantum wire [10,20-23]. Since we will 
use lattice models and the Hartree-Fock (HF) approxima- 
tion for dealing with interactions in our numerical stud- 
ies, we will first discuss those topics in Sec. III. In that 
section, we will also show how the Born approximation 
for scattering can qualitatively explain the dependence 
of the conductance on the wire length which we obtain 
from the RG equations in Sec. II. 

The NEGF formalism will be briefly described in Sec. 
IV. The advantage of this method is that it treats the in- 
finitely extended leads (reservoirs) in an exact way, and 
it can be used for all values of the wire length and the 
impurity strength. However, it is accurate only for weak 
interactions between the electrons because we are forced 
to use a HF approximation for dealing with the interac- 
tions (for reasons which will be explained below) . We will 
use a lattice model for both the wire and its leads. The 
Hamiltonian will have hopping terms in the wire and in 
the leads, and a density-density interaction (on-site for 
spin- 1/2 electrons, and between nearest neighbor sites for 
spinless electrons) only in the wire. 

In Sec. V, we will describe our results for spin-1/2 and 
spinless electrons for different values of the impurity po- 



tential and the interaction strength, and we will compare 
our results with those obtained by the RG analysis. For 
the case of spinless electrons, we find that the agreement 
between the RG and numerical results is excellent if we 
follow a certain procedure. In Sec. VI, we will make 
some concluding remarks. 



II. RENORMALIZATION GROUP EQUATION 
FOR SCATTERING FROM A POINT 

There are several ways of studying the renormalization 
group (RG) evolution of the scattering from one or more 
points in a one-dimensional system of interacting elec- 
trons. One can use the technique of bosonization [7,8,24], 
a fermionic RG method [13-15], and the functional RG 
method [16-19]. Since our numerical calculations use a 
HF approximation, the method of Refs. [13,14] will be 
the most useful for us. Before considering the HF ap- 
proach, however, we will briefly discuss the RG equations 
obtained by bosonization for spinless electrons. 



A. Spinless electrons 

Let us first consider the Hamiltonian for non- 
interacting electrons in the presence of a (5-function im- 
purity placed at the origin, 



H = 



2m 



+ V5{x) 



(2) 



It is easy to check that for plane waves incident either 
from the left or from the right with wavenumber k, the 
reflection and transmission amplitudes are given by [25] 



r(k) 
and t(k) 



imV 



h 2 k 



imV 



H z k 



h k + imV 



(3) 



If the Fermi energy of the electrons is given by E F = 
h 2 kp/(2m) (where kp is the Fermi wavenumber, and 
vp = Tikp/m is the Fermi velocity), then the conduc- 
tance G for spinless electrons at zero temperature is given 
by e 2 /h times \t{kp)\ 2 . For \V\ « hvp, we see that 
\t{k F )\ 2 = 1 - (V/hvp) 2 up to order V 2 . (We will usu- 
ally set Planck's constant h = 1). Eq. (2) has a bound 
state if V < 0, but this does not play any role in the RG 
analysis described below. 

Now let us introduce interactions between the elec- 
trons. We assume a density-density interaction between 
spinless electrons of the form 

Hint = ^ J J dxdy p{x) U(x-y) p(y) , (4) 

where the density p is given in terms of the second- 
quantized electron field ^(x) as p — ^f^. The electron 



2 



field can be written in terms of the right and left moving 
fields and "J" l (whose variations in space are governed 
by wavenumbers much smaller than hp) as 



*(x) = ^ R e 



ikp-x 



-ikpx 



(5) 



If the range of the interaction U (x) is short (of the order 
of A), such as that of a screened Coulomb repulsion, the 
Hamiltonian in (4) can be written in the form 



Hut = .92 J dx *Jj* fl *i,*L , 



(6) 



where gi is related to the Fourier transform of U(x) as 
.92 = U(0) — U(2kp). It is convenient to define the di- 
mcnsionless constant 



a = 



92 

2-jtvf 



(7) 



A system of interacting electrons in one dimension such 
as the one introduced above is described by Tomonaga- 
Luttinger liquid (TLL) theory. The low energy excita- 
tions of a TLL are particle-hole pairs which are bosonic 
in nature and have a linear relation between energy and 
momentum. For spinless electrons, the low energy and 
long distance properties of the TLL are governed by three 
quantities, namely, the velocity v of the low energy exci- 
tations, a dimensionless parameter K which is related to 
the interaction strength, and the Fermi wavenumber Uf- 
For the model described above, we find that [24] 



and K 



v F (1 -a 2 ) 1 / 2 

,1 — 1/2 



(8) 



Thus K = 1 for non- interacting fermions. For weak in- 
teractions, v = vf and K = 1 — a to first order in a. In 
this paper, we will be interested in the case in which the 
interaction is weak, i.e., |a| << 1. 

It turns out that in the presence of interactions, the 
impurity strength V effectively becomes a function of a 
length scale I, and satisfies a RG equation. On bosonizing 
the TLL theory [7,8,24], we obtain the equation 



dV 
dlnl 



(1 - K) V , 



(9) 



to first order in V, i.e., this is valid in the weak barrier 
limit. In the strong barrier limit, which implies \V\ » 
Vf, bosonization leads to a different RG equation 



d(l/V) 
d\nl 



1 



= (!--) (1/V) 



(10) 



which is valid to first order in 1/V. We see that V — 
(V — oo) is a stable fixed point if K > 1 (K < 1 
respectively). Since the above RG equations are not valid 
if \V\/vp is of order 1, we cannot conclude from them 
whether or not there is an intermediate fixed point. 



For a given value of the wire length L and temperature 
T, the RG equations can be used to compute the conduc- 
tance analytically as follows. We begin at a short length 
scale (of order A) with an initial value of V; we then 
use Eq. (9) or (10), depending on whether V is small or 
large, to follow the evolution of V with the length scale 
I. The RG flow stops when I reaches a distance of the or- 
der of L or Lt, whichever is smaller. When the RG flow 
stops, we take the value of V obtained at that point, and 
compute the conductance G in terms of the transmission 
amplitude given in Eq. (3). [We will see later that the 
RG flow of the conductance does not stop abruptly at one 
particular length scale. There is an intermediate range of 
length scales where the conductance continues to evolve 
slowly] 

Let us now discuss the second method for obtaining the 
RG equations, namely, using a HF approximation for the 
case of weak interactions [13-15]. This method directly 
gives an RG equation for the scattering matrix which 
is produced by the impurity. The idea is that reflec- 
tion from the impurity leads to an interference between 
the incoming and outgoing electron waves; this leads to 
Friedel oscillations in the density with an amplitude pro- 
portional to r and a wavelength given by A = ir/kp. 
In the presence of interactions, these density oscillations 
cause the electrons to scatter; these scatterings therefore 
renormalize the scattering caused by the impurity. The 
RG equations obtained by this method are given by 



dt 



a t r ' 



d\nl 

dr l2 
and — — - = a r t 
ami 



(11) 



to first order in the interaction parameter a. (Eq. (11) 
is consistent with the unitarity of the scattering matrix, 
namely, |r| 2 + |t| 2 = 1 and r*t + t*r = 0). The solution 
of Eq. (11) is given by 



W)\ 2 



\t(d)\ 2 {d/lf 



i - W)\ 2 



\t{d)\ 2 (d/l) 



2a 



fl2) 



where d is a short distance scale which is of the order of 
A. However, we will see at the end of Sec. V C that it is 
necessary to choose d to be substantially larger than A in 
order to obtain a good fit between the numerical results 
and the expression in Eq. (12). [At I — ► oo, Eq. (12) has 
fixed points at \t\ = and 1, depending on the sign of a; 
there is no fixed point at intermediate values of t.] 

We can check that Eq. (11) is equivalent to Eqs. (9) 
and (10) if a is small and \t\ is close to 1 and respec- 
tively. If a is small, the Luttinger parameter K ~ 1 — a. 

We note that Eq. (11) is complementary to Eqs. (9) 
and (10) obtained from bosonization in the following 
sense. Eq. (11) is valid for all values of t but only small 
values of a. On the other hand, Eqs. (9) and (10) are 
valid for arbitrary values of the interaction parameter a 
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but only for small values of V or l/V, i.e., only when 
either the reflection amplitude r or the transmission am- 
plitude t is small. 



B. Spin- 1/2 electrons 

We now turn to the more realistic case of electrons 
with spin. For this situation, we will only discuss the 
RG equations obtained from the HF approximation for 
the case of weak interactions [13,14]. We again consider 
a (5-function impurity at one point in the TLL, and a 
density-density interaction of the form given in (4) ; the 
density is now given by 



(13) 



where f and J, denote spin-up and spin-down respectively. 
Introducing right and left moving fields \tfl j(T and $ l,o 
as in Eq. (5), with a =\ or J., we obtain the interaction 
Hamiltonian 

ffint = 

J dx [ 9l *L*L + »2*L*L' *^ 



It is interesting to note that unlike in the spinless case, 
Eq. (16) allows the possibility of \t\ increasing for a while 
and then decreasing (or vice versa) [13]. This is due to 
the flow of the couplings g\ and g 2 ; it happens if U(0) — 
2U(2k F ) and 17(0) - (l/2)U(2k F ) have opposite signs. 
This is precisely the situation for the Hubbard model in 
which U(0) and U(2k F ) are both equal to the Hubbard 
parameter U. 



III. LATTICE MODELS AND THE 
HARTREE-FOCK APPROXIMATION 

Although the system of interest may be defined in the 
continuum, it is convenient to approximate it by a lattice 
model in order to do numerical calculations. (We should 
of course ensure that physical quantities like the wire 
length are much larger than the lattice spacing, so that 
the lattice approximation does not introduce significant 
errors). Let us discuss the form of the lattice Hamil- 
tonians that we will consider. The Hamiltonian has a 
hopping term 



H 



7 S t 



(17) 



(14) 

where g\ — U(2k F ), and g 2 — 34 = U(Q). (We have 
ignored umklapp scattering terms here; they only arise if 
the model is defined on a lattice and we are at half- filling) . 
It is known that g\ , g 2 and g^ satisfy some RG equations 
[26] ; the solutions of the lowest order RG equations are 
given by [13] 



9i (0 



U(2k F ) 



U(2k F ) 

17V F 



92(1) = U(0) -\ U(2k F ) + i 



9i(l) = U(0) . 



U(2k F ) 



U(2k f 



■ In I 



(15) 



Next, we can use the HF approximation and the exis- 
tence of Friedel oscillations in the density to derive the 
RG equation of the transmission amplitude t. We find 
that [13,14] 



dt 
dhd 



92(1) - 2gi(0 
2ttv f 



t\rf 



(16) 



Due to the RG flow of the couplings g\ and g 2 , the solu- 
tion of Eq. (16) is more complicated than the expression 
for spinless fermions given in Eq. (12) [13]. We will 
therefore not attempt to make a quantitative compari- 
son between our numerical results and the RG equations 
for the case of spin-1/2 electrons. 



where c n _ a annihilates an electron with spin a at site n, 
and the hopping amplitude 7 will be taken to be positive 
for convenience. Next, we will place an impurity at one 
site, say, n = 0, so that 



(18) 



Finally, we will introduce an interaction between the elec- 
trons in the wire (but not in the leads which will be con- 
sidered in the next section). For the spin-1/2 case, the 
interaction will be taken to be of the Hubbard form, 



Hint = U C 1,T C «,T C l| C «< 



(19) 



For the case of spinless electrons, the spin index a will be 
dropped in Eqs. (17) and (18), and the interaction will 
be taken to be between nearest neighbor sites, 



IS 



T +l<-n+l 



(20) 



(For spinless electrons, we cannot have an on-site inter- 
action since (c) n c n ) 2 = c„). 

In the absence of the impurity and interactions, the 
energy is related to the wavenumber as 



E(k) 



27 cos A; 



(21) 



where — 7r < k < it. (We have set the lattice spacing 
equal to 1). The velocity is given by v = dE/dk = 
2jsink. If the chemical potential is given by fi, with 
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—27 < p < 27, the ground state is the one in which 
the states with momenta — kp to kp arc filled, where 
/1 = — 2-f cos kp, with < kp < n. It is convenient 
to define the Fermi energy as the difference between the 
chemical potential and the bottom of the band, namely, 
E F = p + 2-f = 2j(l - coskp). 

Let us now consider the effect of the impurity placed 
at n — (we are still ignoring the interactions). The 
reflection and transmission amplitudes are given by 



r(jfe) 
and t(k) 



iV 



27sinfc - 
27 sin k 



iV 



27 sin k 



iV 



(22) 



In addition, there is a bound state if V < 0; the bound 
state energy is given by E\, = —2-fcoshkb, where V = 
— 27sinhfcb, with fc& > 0. The normalized bound state 
wave function is given by 



tpb,n = (tanhfc fc ) ' e 



1/2 -k b \n\ 



(23) 



for all values of n. 

It is interesting to compute the electron density p n = 
Pn,] = Pn,i, as a function of the site index n. The con- 
tribution from the scattering states is given by 



o n = I g [|e- ife H + r(k)e^\ 2 + \t(k)\ 2 

"' [ r(k)e l2k ^ + r*(k)e- l2k ^ } , 

(24) 



= P 



\e 

kF dk 
2?r 



where p = kp/n is the value of the density very far from 
the impurity. In the limit \n\ — > 00, we find that the 
density has oscillations given by 



Pn ~ P 



47r|n| 



r(kp)e 



i2k F \n\ _ r *(/ CF ) e - i2fc Fl™ 



(25) 



Now we will introduce interactions and consider the 
HF approximation for the cases of spin- 1/2 and spinless 
electrons. This will give us another way of understanding 
the RG equations discussed in Sec. II. 



A. Spin-1/2 electrons 

For the interaction given in Eq. (19), the HF approxi- 
mation takes the form 

ffint.HF = U [Pn (4,T C «,t + C n,i C "J ~ Pn } , 

n 

(26) 

where we have set < ^c n ^ >—< <y n ^Cn,i >= Pn- 
[This is called the restricted HF. One can also make 



an unrestricted HF approximation in which one allows 
< c„ iT c„ ;T > to differ from < c^c^ >.] 

At this point, we would like to impose the condition 
that the interaction should have no effect on the conduc- 
tance if there is no impurity, i.e., if p n is equal to the 
constant p at all sites. The reason for imposing this con- 
dition is that we do not want the interactions by them- 
selves to lead to any scattering; that would change the 
conductance from the ideal value of 2e 2 /h in the absence 
of any impurities which is an undesirable effect. (This 
will become clearer in the next section when we describe 
our model for the wires and the leads). We will therefore 
modify Eq. (26) to the form 

ffint,HF = U ^ (p n - p) (4 ;T C„, T + cl d C nd ) , 
n 

(27) 

where we have dropped the constant p\ which turns out 
to have no effect on the conductance. (The terms pro- 
portional to p in Eq. (27) are equivalent to adding a 
chemical potential). 

To gain an insight into the effect of Eq. (27), let us 
consider the Born approximation. In order to use this 
approximation, we will assume here that both V and U 
are small. The total potential seen by either spin-up or 
spin-down electrons at site n is 



V n = VS nfi + U{ Pn - p) 



(28) 



For small values of |V|/i|f (where vp = 2-f sin kp is the 
Fermi velocity), r(kp) = —iV/vp from Eq. (22), and we 
have 



V n = V8 n n 



vu 



A-KVp\n\ 



( e i2fc F |n| + e -i2fc F M) j (29) 



where the second term is valid only for \n\ » ir/kp. 

The Born approximation for the reflection amplitude 
in one dimension for a lattice model is given by 



r B (k F ) - -— VV, 

VF L — ' 



(30) 



If the region in which the electrons interact with each 
other only extends from n=— L/2 to L/2, then for L >> 
n/kF, we get 



TB{kp) = - 



iV 



VF 



[1 



-i^ln^ 

2-kvf 2 



(31) 



If U > 0, we see that the conductance G = (2e 2 /h)[l — 
\rB{kp)\ 2 ] increases with L till it reaches a maximum 
when ln(L/2) ~ 2ttvf/U. However, the Born approxima- 
tion that we are using here cannot be trusted up to such 
large length scales because we are not self-consistently 
modifying the densities at the different sites in response 
to the renormalization of the reflection amplitudes (pro- 
duced by the interaction U). Because of this lack of self- 
consistency, we can only trust the Born approximation 
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results up to length scales where the reflection amplitude 
has not changed very much from the value of —iV/vf 
that it has in the non-interacting theory. 

We note that Eq. (31) is consistent with Eq. (16) 
for small values of V and (U/2ttvf)^1, since (32(0 — 
2gi(i))/2irvp = — U/2nvp for such values of I. 



Hence, 



47r|n| 



- r*{k F )e 



-i2k.F \n\ 



(37) 



B. Spinless electrons 

For the interaction given in Eq. (20), the HF approxi- 
mation takes the form 

U \ , 4. + 

«int,HF = 2-^i I P n+1 n ° n + PnC n+1 C n+ i 

n 

— < c|jC„+l > c n+1 c n 

- < C l+l C n > cUn+l } , (32) 

where p n =< $ n c n >, and we have dropped some con- 
stants. In the absence of any impurities, p n — kp/ir and 
< c\ l+1 c n >= (sinfci?)/7r for all values of n. 

Once again, we impose the condition that the interac- 
tion should have no effect on the conductance if there is 
no impurity. We therefore modify Eq. (32) to the form 

ffint.HF = y J] [{Pn+1 + Pn-1 ~ 2p) C n C n 
n 

, t sinfc F t 

- (< c\Cn+l > ) c! n+1 C n 

7T 

, + sinfcF, + 1 

- (< C n+l C n > ) c\ Cn + X } , 

7T 

(33) 

where p = kp/ir as before. 

Let us again use the Born approximation to understand 
the effect of Eq. (33), assuming that both V and U are 
small. If one adds a perturbation of the form 

#pcrt = t V " C " C ™ 

n 

+ Vn+l/2 C f n+1 C n + V* +1/2 clcn+l } 

(34) 

to the Hamiltonian in Eq. (17), the Born approximation 
for the reflection amplitude is given by 



r B (k F ) = -— V [V„e f 
vf — 



2kp-n 



+(V„+i/2 + K +1/2 ) e 



:fc F (2n+l)i 



(35) 



In our case, 



U 



Vn = VS n fi + ^" I Pn+1 + Pn-1 - 

U , + sinfc F 

Vn+l/2 = - "2 I < c n c n+l > ~ 



2p 



where the second term is only valid for \n\ — > 00, and 

iU 



K+1/2 



87r|n| 



[r(kp) e i2k P\ n \+isgn(n)kF 
— T*{kp) e - i2k FW\-iagn{n)k F 



(38) 



for |n| — > 00; here sgn(n) = n/|n|. For small values of 
jl^l/^F, f(kp) = —iV/vp. Substituting everything in the 
Born approximation (35), and assuming that the region 
in which the electrons interact with each other extends 
from n = —L/2 to L/2, where L » ir/kp, we find that 



iV 

r B (k F ) = - — [ 1 
vf 



JL-(1 -oos(2* F )) In I 



(39) 



(36) 



Eq. (39) is consistent with Eq. (11) for small values of 
V and (U/2irv F )lnl, since a = (U/2irv F )[l - cos(2fc F )]. 



IV. NON-EQUILIBRIUM GREEN'S FUNCTION 
FORMALISM 

In this section, we will introduce the NEGF formalism 
which will allow us to study the conductance of a wire 
with any length, both short (of the order of ir/kp) and 
long (where a continuum description and RG analysis 
may be expected to be reliable). 



A. Self-energy, density and conductance 

An important concept in the NEGF formalism is a 
"self-energy" which describes the amplitude for an elec- 
tron to leave or enter the wire from the leads (reser- 
voirs) which are maintained at some chemical potentials 
and temperature [10,20-23]. The self-energy is a non- 
Hcrmitian term in the single-particle Hamiltonian of the 
wire. To see how this arises, let us begin by modeling 
one of the reservoirs by a tight-binding Hamiltonian 

H = -7 E ( C n,* c n+l,a + 4+1,<t c ™^ ) • (40) 
n.cr 

The energy of an electron in the reservoir is related to 
its wavenumber as E = —2"fcosk. The reservoir has 
a chemical potential p and an inverse temperature () — 
l/(fcsT). The reservoir is semi- infinite; the last site at 
one end of the reservoir is coupled to the first site of the 
wire by a hopping amplitude 7' (see Fig. 1). 
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FIG. 1. (Color online) Picture of a semi-infinite reservoir 
going from site number to -co and a wire beginning from 
site number 1. 



The reservoir gives rise to a self-energy at the first site 
of the wire of the form [10,22] 



where 



£(£) = a(E) ]T 4 )CTCl „ 



a(E) = - — e lk 

7 



i 2 - k ■ 

— — e — it] 

7 



for -27 < E = -27 cos fc < 27 (with < k < tt), 

<t(25) = 

for E = -27 cosh fc < -27 (with fc > 0), and 
<t(£7) = J— e ~ fc - i?7 



(41) 



(42) 



(43) 



(44) 



for E — 27 cosh fc > 27 (with fc > 0). Here r\ is an 
infinitesimal positive number which appears only if E 
lies outside the range [—27,27]. For —27 < E < 2j, the 
self-energy already has an imaginary piece, so it is not 
necessary to add an in term. 



Left 
Reservoir 

ML 



7 71 71 71 
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FIG. 2. (Color online) A quantum wire system with a 
iV-site wire in the middle and semi-infinite reservoirs on its 
left and right. 

Now let us consider the complete system which consists 
of a wire with N sites in the middle, and reservoirs on 
its left and right as shown in Fig. 2. If the wire is also 
modeled by a tight-binding Hamiltonian (with a hopping 
amplitude 71 ) , the Green's function of the wire at energy 
E is given by a N x N matrix 

G{E) = [EI - H - E L (E) - Z R (E) p 1 , 





+ 


X R (E) 


= 


rr(E) 









7i 





7i 








7i 





7i 









7i 



(45) 



7i o{E) j 



(The effects of the reservoirs on the wire is completely 
taken into account by the self-energy terms). We would 
like to emphasize that the relation between E and fc is 
given entirely by the dispersion in the reservoirs (we as- 
sume that the dispersion is the same in both the reser- 
voirs), and not on the form of the Hamiltonian H inside 
the wire. 

The density matrices due to electrons coming in from 
the left and right reservoirs are given by [10,20-22] 



/ S GrLGt/L and / S Grv?t h 

respectively, where T a (E) = i(E a — £t), 



fa(E) 



(46) 



(47) 



is the Fermi function, and \i a is the chemical potential 
in reservoir a. We note that if E lies outside the range 
[—27,27], the matrix T a (E) ~ 2rj is infinitesimal. In 
that case, the density matrix can still receive a contribu- 
tion from certain states; these are typically bound states 
which have a discrete set of energies £&. The reason that 
such states can make a finite contribution in Eq. (46) 
even though r a is infinitesimal is that G takes the form 
1/(E -E b + irj), and 



Lim 



277 



= 2 tt 5{E - E b ) . (48) 



{E-E b f 

Finally, the current is given by the expression 

j dE trace(Gr L G t r fl ) [f L (E) - f R (E)] , 

(49) 



e 
h 



and the conductance isG = e — Mr) ■ If E 

lies outside the range [—27,27], Tl(E) and Tr(E) are 
both infinitesimal, and the current does not get any con- 
tribution from states lying in that energy range. (The 
difference between the density matrix and the current 
for such states is that the density matrix has only one 
factor of 77 in the numerator while the current has two 
factors of 7/). 

We will be interested below in the case of linear re- 
sponse, i.e, the limit \il — * fJ-R = M- In that limit, the 
Fermi wavenumber kp is given by fi = — 2^/coskF- Fur- 
ther, the conductance takes the form 
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P 2 z-27 a 

G = — / tracefGTLG^^) -f . 

(50) 

(We have to multiply the above expression by 2 for spin- 
half electrons). 

B. Interactions 

Let us now consider how interactions can be studied 
within the NEGF formalism. Note that this formalism 
works with a one-particle Hamiltonian, e.g., the self- 
energy in Eq. (41) is given in terms of the energy of 
a single electron which is entering or leaving the wire. 
Hence, the only way to deal with interactions is to do 
a HF decomposition as shown in Eqs. (26) and (32); 
these take the form of corrections to the on-site chemical 
potential or the hopping amplitude. 

However, there is a difficulty in using the expressions 
in (26) or (32). We want to have interactions between the 
electrons only in the wire, not in the reservoirs. The HF 
decompositions shown in (26) and (32) will then lead to 
one-body terms only in the wire; these terms will back- 
scatter electrons coming from the reservoirs, and hence 
reduce the conductance from its ideal value of 2e 2 /h or 
e 2 /h. We would like to ensure that there is no such scat- 
tering if there are no impurities inside the wire and if the 
leads connect adiabatically to the wires, i.e., if the entire 
system with wire and reservoirs is translation invariant. 
We know that two-body scatterings between electrons 
conserve momentum and therefore do not affect the con- 
ductance in the absence of impurities and scattering from 
the lead-wire junctions. We therefore modify the form of 
the HF decomposition from Eqs. (26) and (32) to Eqs. 
(27) and (33) respectively; the extra terms in those equa- 
tions ensure that the HF terms vanish identically if there 
are no impurities. [An alternative way of ensuring adia- 
baticity between the leads and the wires is to turn on the 
interaction strength U smoothly from zero in the leads 
to a finite value inside the wire. This, however, is harder 
to implement numerically. Subtracting the mean values 
as in (27) and (33) is easier to implement, and it serves 
the same purpose because the subtracted quantities ap- 
proach zero near the ends of the wires (this will become 
clear from the numerical results presented below).] 

In the presence of interactions, we will implement a 
self-consistent NEGF calculation as follows: 

• Start with the Hamiltonian with no interactions, 
and calculate the density matrix. The diagonal el- 
ements of the density matrix give the densities at 
different sites. 

• Use the HF approximation to compute the Hamil- 
tonian with interactions, and use that to calculate 
the density matrix again. 



• Repeat the previous step till the density changes 
no further, i.e., between two successive iterations, 
the maximum change in density on any lattice site 
is less than about 10~ 4 . 

• Use the converged density to compute the conduc- 
tance. 

V. NUMERICAL RESULTS 

In this section, we will present the results obtained by 
the NEGF formalism, first for non-interacting electrons, 
and then for interacting spin- 1/2 and spinless electrons, 
using the procedure described in Sec. IV B. We will study 
the dependence of the conductance G on the wire length 
L, the inverse temperature 0, the impurity strength V, 
and the interaction parameter U. Finally, we will make 
a quantitative comparison between our numerical results 
and the RG equations for the case of spinless electrons. 

We present some details about our choice of parameters 
for the calculations. We always take the wire to have an 
odd number of sites, and the impurity to be on the middle 
site. The hopping amplitudes in the reservoirs and in 
the wire is chosen to be the same (7 = 71 = 7' = 1); 
the values of the inverse temperature = l/(fcsT) is 
quoted in units of I/7. (We set both Ti and the lattice 
spacing equal to 1). We choose fc F = n/10; hence, vp = 
2sinfc F = 0.618, fj, = -2cosfc F = -1.902, and E F = 
2(1 — cosfc F ) = 0.098. The thermal coherence length is 
given by Lt = 0.618/3. In the energy integrations for the 
current and conductance, we take the energy step size to 
be dE = 2 x 10~ 4 . In the integration outside the energy 
range [—27, 27], we take the quantity 77 to be 5 times dE. 
Finally, in all the figures, the conductance G is expressed 
in units of 2e 2 /h and e 2 /h for the spin-1/2 and spinless 
cases respectively. 

Our choice of parameters was dictated partly by ex- 
periments on quantum wires in semiconductor hetero- 
junctions [1-6], and partly by our numerical limitations 
(which prevent us from going to very large wire lengths) . 
Experimentally, the ratio of the wire length to the de 
Broglie wavelength of the electrons (Lfc F /7r) ranges from 
about 20 to 200, the ratio of the wire length to the inverse 
temperature L/0 ranges from about 1/2 to 10, and the 
parameter U/(2ttvf) appearing in the Eq. (31) is about 
0.2 - 0.3 [9]. In our calculations, Lkpf-K goes from about 
1 to 30, L/0 goes from about 1/40 to 6, and the interac- 
tion parameters take the values U / (2ttvf) — 0.0773 in the 
spin-1/2 case (Eq. (31)) and (U/2nv F )[l - cos(2/c F )] = 
0.0148 in the spinless case (Eq. (39)) for U = 0.3. The 
much smaller value of the interaction parameter in the 
spinless case leads to smaller changes in the conductance 
in that case as we will see. 
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A. Non-interacting electrons 



Let us first discuss some properties of a system in which 
there is a (5-function impurity of strength V in the mid- 
dle of the wire, and there are no interactions between 
the electrons. Fig. 3 shows the Friedel oscillations in the 
density < c}jC„ > for V = 0.3 for two different tempera- 
tures. Note that the density at the middle site is much 
lower than the mean value of p = kp/n — 1/10 due to 
the repulsive nature of the impurity. We also observe 
that the oscillations die out beyond a length scale of or- 
der Lt = 0.618/3 in the lower figure. This is a simple 
illustration of the idea of a thermal coherence length. 



(3=300 
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0.09 - 
Q_ 
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FIG. 3. (Color online) Friedel oscillations in the density of 
non-interacting electrons caused by an impurity with V = 0.3 
placed in the middle of a wire with 151 sites. The upper and 
lower figures have (3 = 300 and 50 respectively. 

Fig. 4 shows the Friedel oscillations in the density 
for V = -0.3 for (3 = 300. The density in the middle 
site is now much higher than the mean value; this is 



(3=300 




FIG. 4. (Color online) Friedel oscillations in the den- 
sity of non-interacting electrons caused by an impurity with 

V = —0.3 placed in the middle of a wire with 151 sites and 
/3 = 300. 

because of the presence of a bound state whose wave 
function has a peak at that site. Indeed, in the absence 
of interactions, the difference in the site densities for V — 
0.3 and —0.3 is given entirely by the bound state. We 
see from Figs. 3 and 4 that the density difference at the 
site of the impurity is about 0.15; this agrees well with 
the value of tanhfc & = 0.148 given by Eq. (23). (For 

V = -0.3, we have k b = 0.149). 

Finally, Fig. 5 shows the quantity p\ =< c3 n c n +\ > for 

V = 0.3 and —0.3 for (3 — 300. There is a strong similar- 
ity between these and the density plots shown in Figs. 3 
(upper figure) and 4. This is because the electrons at the 
Fermi energy (which dominate the long distance proper- 
ties of the system) have a wavelength which isir/kp = 10 
times longer than the lattice spacing; hence there is very 
little difference between < cJjC„ > and < cJjC„+i >. 

For non-interacting electrons, the conductance de- 
pends on the temperature, but not on the length of the 
wire. For a 5-function impurity of strength V, the con- 
ductance for spinless electrons is given by 



G = 



dE T{E) 



T{E) = 



1-2 

4 sin 2 k 
4 sin 2 k + V 2 



[2 cosh§(£; 



■A*)] 5 



(51) 



where E = -2 cos A;. For k B T « E F = fi + 2, Eq. (51) 
has a Sommerfeld expansion of the form [27] 



G 



T(p) 



7T 



(k B Tf T"(jx) 



(52) 



For k F = tt/10 and V = 0.3, T{p) = 0.809 and 
(it 2 /6)T"(p) — 42.1; this implies a significant temper- 
ature dependence of the conductance. Fig. 6 shows that 
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FIG. 5. (Color online) Plots of pi =< cj,c n +i > for 
non-interacting electrons in the presence of an impurity with 
V — 0.3 (upper figure) and —0.3 (lower figure) placed in the 
middle of a wire with 151 sites, with /3 = 300. 



the conductance changes appreciably with (3 till (3 reaches 
about 100. We have to keep this in mind when study- 
ing the temperature dependence of the conductance of a 
system of interacting electrons. 



B. Spin-1/2 electrons 

We now consider the spin-1/2 model with an on-site 
interaction between the electrons. In the absence of in- 
teractions, the conductance at zero temperature is inde- 
pendent of the wire length and is given by Eq. (52), with 
a factor of 2 for spin; namely, 



4 sin ki 



c - 2e2 

h 4 sin kp 



V 2 



(53) 



For k F = tt/10 and \V\ = 0.03, G = {2e 2 /h) 0.809. 



FIG. 6. (Color online) Conductance versus /3 for 
non-interacting electrons in the presence of an impurity with 
V = 0.3 placed in the middle of a wire with 201 sites. 



In Fig. 7, we show the conductance as a function of 
the wire length L for three different value of P, with 
U = 0.3 in the upper figure and U = —0.3 in the lower 
figure. The trends in Fig. 7 are in accordance with Eq. 
(31). Firstly, the conductance increases with the length 
scale if U > and decreases ii U < 0. For U = 0.3, 
we have U/(2irvp) — 0.0773; we expect the HF approxi- 
mation to be reasonable for such a small value. [Due to 
the RG flows of gi and gi in Eq. (16), the conductance 
should start decreasing for positive U and increasing for 
negative U at a very large length scale of the order of 
L ~ exp(l/0.0773) ~ 420000, provided that L T is also 
larger than 420000. Such a wire length is beyond our 
numerical capability.] 

Secondly, for fixed f3, G increases till L reaches a value 
of the order of Lt = 0.618/? beyond which G stops chang- 
ing. This happens because there is complete thermal dc- 
coherence once L exceeds Lt; hence G does not vary 
any more with L in that regime. Thirdly, the conduc- 
tance for two different values of /3, say, /3i and /?2 (where 
Pi < P2), start separating from each other at a value 
of the wire length L c such that L c / Pi is much less than 
0.618. This happens because if the inverse temperature 
is Pi or higher, there is complete thermal coherence, and 
G does not depend on P any more. Thus, for any temper- 
ature P, there is a range of wire lengths going from L c up 
to Lt where there is partial thermal coherence, and the 
conductance evolves slowly in this range of lengths. Wc 
therefore conclude that the RG flow does not suddenly 
stop at a length scale which is the smaller of the wire 
length L and the thermal coherence length Lt] rather, 
the flow continues to occur (but slows down) in an inter- 
mediate range of length scales. 

Finally, Fig. 7 shows that even for the smallest values 
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FIG. 7. (Color online) Conductance versus wire length for 
interacting spin-1/2 electrons for three different values of (3, 
for V = 0.3. U = 0.3 and -0.3 in the upper and lower figures 
respectively. 



of L, the conductance differs from the non- interacting 
value of 0.809 by an appreciable amount. This is due to 
the large deviations of the density from the mean value at 
the sites close to the impurity. In that region, the density 
does not satisfy the l/\n\ form given in Eq. (25). Since 
that form is intimately tied to the RG results (the sum 
over l/\n\ gives ln(L)), we do not expect the numerically 
obtained value of the conductance to agree well with the 
RG analysis for short wire lengths. 

In Fig. 8, we show the conductance as a function of (3 
for three different value of L, with U = 0.3. Once again, 
we see the same trends as in Fig. 7. Namely, for any 
two values of L, there is a (3 where the values of G start 
separating from each other. Then there is a higher value 
of (3 where G stops changing. The region between the 
two is where there is only partial thermal coherence, and 




FIG. 8. (Color online) Conductance versus (3 for interact- 
ing spin-1/2 electrons for three different wire lengths, for 

V = 0.3 and U = 0.3. 

G varies relatively slowly in that region. 

In Fig. 9, we show the conductance as a function of 

V for three different values of U, with L = 101 and j3 — 
100. Note that for U = 0, the conductance is an even 
function of V given by Eq. (53). But in the presence 
of interactions, G is no longer precisely an even function 
of V , particularly for large values of V . This is more 
clearly visible in Fig. 10 which shows the conductance as 
a function of U for V = 1 and -1, with L = 101 and (3 = 
100. For any value of U, we find that the conductance 
deviates less from the non-interacting value for V > 
compared to V < 0, for the same value of |V|. For small 
values of U, this can be explained as follows. We saw in 
Figs. 3 and 4 that the density at the site of the impurity 
deviates from the mean density of p — 1/10 by a smaller 
amount for V > compared to V < 0; the numerical 
values for po — 1/10 are given by -0.046 and 0.095 for 

V = 0.3 and -0.3 respectively. The effective impurity 
strength is given by 



V + U (p - p) 



(54) 



which equals 0.3 - 0.046C7 and -0.3 + 0.095J7 for V = 
0.3 and —0.3 respectively. For U small and positive, 
this means that the magnitude of the effective impu- 
rity strength and therefore the reflection probability is 
larger for V = 0.3; hence the conductance is smaller for 

V = 0.3. The opposite statement is true if U is small and 
negative. To conclude, the different values of the density 
deviation at the impurity site for positive and negative 
values of V are responsible for the asymmetry in G as a 
function of V. (This asymmetry is more prominent for 
the case of spinless electrons as we will see below). 

The fact that the deviation of the conductance from 
the non-interacting value is smaller for V > than for 

V < is clearly a result of the bound state which raises 
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FIG. 9. (Color online) Conductance G versus V for inter- 
acting spin-1/2 electrons for U = 0.3, and -0.3, with L = 101 
and = 100. 

the density at the impurity site for V < 0; hence this is a 
short distance effect. [The contribution of the scattering 
states to the reflection probability |r| 2 is an even function 
of V as one can see from Eq. (31). The RG equations 
which are derived from continuum theories take into ac- 
count only long distance effects; these come from the 
scattering states only.] 

In Fig. 11, we show the Friedel oscillations for three 
different values of U with V = 0.3, L = 101 and = 100. 
We see that a repulsive interaction (U > 0) suppresses 
the Friedel oscillations, while an attractive interaction 
(U < 0) with the same magnitude enhances the oscilla- 
tions by a much larger amount. This is because attrac- 
tive interactions tend to lead to the formation of a charge 
density wave; if a charge density wave is already present 
(due to the impurity), attractive interactions enhance it. 
Since the scattering from the Friedel oscillations renor- 
malize the scattering from the impurity, this difference in 
the magnitude of the oscillations helps to explain why the 
deviation of the conductance G from the non-interacting 
value is larger for U < compared to U > 0, as can be 
seen in Fig. 7. 

For spin-1/2 electrons, we can do a self-consistent HF 
calculation in two ways, restricted and unrestricted. In a 
restricted HF calculation, the site densities of spin-up and 
spin-down electrons are taken to be equal at all stages. 
In a unrestricted HF calculation, spin-up and spin-down 
electrons are allowed to have different densities. We have 
done both kinds of calculations. For the range of the in- 
teraction —0.5 < U < 0.5, we find that they give the 
same results; the spin-up and spin-down densities con- 
verge to the same values even if one begins with different 
initial values for them. Thus we do not find any spin 



FIG. 10. (Color online) Conductance G versus U for inter- 
acting spin-1/2 electrons for V = 1 and -1, with L = 101 and 
= 100. 

density wave within our range of parameters. 

C. Spinless electrons 

We now consider a model of spinless electrons with 
nearest neighbor interactions between the electrons. 
Comparing the forms of Eqs. (31) and (39), we see 
that the behaviors of the spin-1/2 model for U > 
and U < should be similar to that of the spinless 
model for U < and U > respectively. In Fig. 
12, we show the conductance as a function of the wire 
length L for three different value of 0, with V = 0.3 
and U = 0.3. The trends are in agreement with Eq. 
(31), with a = U(l- cos(2fc F ))/(27ri> F ) = 0.0148. For 
fixed 0, G decreases till L reaches a value of the order of 
L T = 0.618/3. 

In Fig. 13, we show the conductance as a function of 
for three different value of U, with V = 0.3 and L = 201. 
As we saw earlier in Fig. 6, G has a significant depen- 
dence on for small values of even for non-interacting 
electrons. It is only when we go to large values of 
that we see the trend expected from the RG equations 
for interacting electrons; namely, as increases, the con- 
ductance decreases if U > and increases if U < 0. 

Since G has an appreciable dependence on even for 
non-interacting electrons, it is easier to consider the de- 
pendence of G on the wire length L at zero temperature 
in order to see how well the numerical results compare 
with the RG expression given in Eq. (12). That equation 
has two parameters, namely, the interaction parameter 
a and the short distance scale d (with its corresponding 
transmission probability |t(d)| 2 ). To begin, let us choose 
a = 0.0148 as given by the analytical expression in Eq. 
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FIG. 11. (Color online) Friedel oscillations for spin-1/2 
electrons for U = 0.5, and -0.5 (upper, middle and lower 
figures respectively) for V = 0.3, L — 101 and (3 = 100. 



FIG. 12. (Color online) Conductance versus wire length for 
interacting spinless electrons for three different values of ft, for 
V = 0.3 and U = 0.3. 



U = - 0.5 




FIG. 13. (Color online) Conductance versus /3 for interact- 
ing spinless electrons for three different values of U, for a wire 
with 201 sites and V = 0.3. 



(39) in terms of U and vf- In Fig- 14, we show a compar- 
ison between the expression in (12) for two values of d, 
namely, d = 11 with \t(d)\ 2 = 0.805 (the dash dot line A), 
and d = 51 with \t(d)\ 2 = 0.795 (the dashed line B), with 
our numerical results shown by astersisks. (The values of 
|£(d)| 2 for d = 11 and 51 have themselves been obtained 
by our numerical calculations). We see that the RG ex- 
pression with d = 51 fits the numerical results better than 
the expression with d = 11. This difference between the 
lines A and B shows that the interactions change the con- 
ductance by an appreciable amount between d = 11 and 
d = 51, and this change is not captured accurately by the 
expression in Eq. (12). [The conductance at short dis- 
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FIG. 14. (Color online) Conductance G versus L for inter- 
acting spinless electrons compared with the RG expression 
in Eq. (12), for V = 0.3, U = 0.3, and (5 = oo. Line 
A (dash dot line) corresponds to d — 11, i(d)| 2 = 0.805 
and a = 0.0148. Line B (dashed line) corresponds to 
d = 51, \t{d)\ 2 = 0.795 and a = 0.0148. Line C (solid line) 
corresponds to d = 51, |t(d)| 2 = 0.795 and a = 0.016. The 
asterisks show the numerical results. 
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FIG. 16. (Color online) Conductance G versus U for inter- 
acting spinless electrons for V = 1 and -1, with L = 101 and 
/3 = 100. 



that even line B corresponding to d = 51 starts deviat- 
ing from the numerical results at large length scales. We 
have therefore tried varying the parameter a also, keep- 
ing d and |£(g?)| 2 fixed at 51 and 0.795 respectively. We 
find that a — 0.016 (the solid line C in Fig. 14) gives an 
excellent fit to the numerical results. 

In Fig. 15, we show the conductance as a function of V 
for three different values of U, with L = 101 and (3 = 100. 
We see that G is not an even function of V in the presence 
of interactions. In Fig. 16, we show the conductance as 
a function of U for V = 1 and -1, with L = 101 and 
(3 = 100. For any value of U, we again find that the 
conductance deviates less from the non-interacting value 
for V = 1 compared to V = — 1. The explanation for 
this is similar to that for the similar phenomenon in the 
spin- 1/2 model, except that the roles of U and —U are 
interchanged in the two models. The larger change in the 
conductance for V — — 1 is again due to the presence of 
a bound state. 



FIG. 15. (Color online) Conductance G versus V for inter- 
acting spinless electrons for U = 0.3, and -0.3, with L — 101 
and P = 100. 

tances is affected substantially by quantities such as the 
large value of p — p at the site of the impurity. The RG 
analysis, on which (12) is based, does not take such effects 
into account.] Since the wavelength of the Friedel oscilla- 
tions is given by A = 10, we conclude that the RG results 
agree reasonably well with the numerical results only if 
we begin the integration of the RG flows from a distance 
which is significantly larger than A. However, we observe 



VI. DISCUSSION 

We have used the NEGF formalism to study the de- 
pendence of the conductance of a quantum wire with 
both spin-1/2 and spinless electrons on various param- 
eters such as the wire length, temperature, impurity po- 
tential, and the strength of the interactions between the 
electrons. The advantage of the NEGF formalism is that 
it can be used to compute the density and conductance 
at any length scale. At large length scales, our numerical 
results agree with those obtained by an RG analysis of 
continuum theories. We find that the trends of the RG 
results can be understood using the Born approximation 
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for scattering from a weak impurity and from the density 
oscillations produced by the impurity. 

Our numerical results differ in detail in two ways 
from those obtained analytically from the RG equations. 
Firstly, the dependence of G on the wire length and tem- 
perature fits the expression in Eq. (12) only if we take 
the starting point of the RG equation to be significantly 
larger than the short distance scale A, and we allow a to 
be a little different from its analytically obtained value. 
Secondly, the conductance is not an even function of the 
impurity potential V (as one expects from the RG anal- 
ysis) ; this is due to the existence of a bound state for an 
attractive impurity. These differences between the nu- 
merical results and the results based on the RG analysis 
seem to be due to effects at short distances, where the 
density deviates significantly from the mean density. 

Before ending, we would like to briefly compare our 
work with that of Refs. [16-18]. Using the functional RG 
technique, the authors of those papers have shown that 
there is excellent agreement at very large length scales 
between the numerical results and the asymptotic scaling 
forms given by the RG equations derived from continuum 
theories. We have not been able to go up to such large 
length scales using the NEGF formalism. However, even 
at the length scales studied by us, our numerical results 
agree quite well with the RG equations if we start at a 
short distance scale which is about 5 times larger than A. 

As mentioned towards the beginning of Sec. V, the 
length scales L and Lt that we have considered are com- 
parable to those studied experimentally. Our observa- 
tions about the short distance effects may therefore have 
implications for fitting experimental data to expressions 
obtained by an RG analysis. 
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